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We study quark and electron EDMs generated by Yukawa couplings in supersym- 
metric models with different gauge groups, using the EDM properties under flavour 
transformations. In the MSSM (or if soft terms are mediated below the unification 
scale) the one loop contributions to the neutron EDM are smaller than in previous 
computations based on numerical methods, although increasing as tan 3 /?. A neu- 
tron EDM close to the experimental limits can be generated in SU(5), if tan/3 is 
large, through the u-quark EDM d u , proportional to tan 4 /3. This effect has to be 
taken into account also in SO(10) with large tan/3, where d u is comparable to the 
d quark EDM, proportional to tan/3. 



1 Introduction 

In supersymmetric theories there are various possibilities of generating an intrinsic electric 
dipole moment (EDM) in microscopic systems. 

First of all, the neutron has already a non vanishing EDM in the Standard Model (SM), due 
to the uneliminable phase in the CKM matrix, or, in other words, to the misalignment between 
the left eigenstates of the two Yukawa matrices, Au and A D , and to the charged current gauge 
interactions between the left-handed quarks, that make the misalignment physically significant. 

If further interactions are present, further misalignments and phases can become significant. 
For example, in the unified extensions of the SM the new gauge interactions make significant 
the misalignments among particles unified in the same gauge multiplet. Nevertheless, without 
supersymmetry, due to the decoupling of heavy particles, these effects are suppressed at low 
energy by powers of the unification mass. 

In the supersymmetric extensions of the SM, one more possible source of flavour and CP 
violation is associated to the possibility that the scalar partners of quarks and leptons are not 
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degenerate and point in different directions in the flavour space compared to the corresponding 
fermions. On the other hand, models are often considered in which the soft supersymmetry 
breaking terms are flavour universal and real at some energy scale Mo. This can happen in 
models in which supersymmetry breaking is communicated to the observable sector by grav- 
ity [1, 2] (in this case M = Mp = reduced Planck mass) or when it is communicated by gauge 
interactions at relatively low energy [3] . These hypotheses allow to keep under control flavour 
changing neutral current processes and EDMs due to the mixings in the fermion-scalar interac- 
tions. In this case, effects coming from the non universal radiative corrections to soft breaking 
parameters have to be ascribed to the Yukawa couplings. 

As said, when the gauge group is unified, some new phases and flavour mixings otherwise 
inexistent or unphysical can become significant. Unlike what happens without supersymmetry, 
if its breaking is communicated to the observable sector at a scale where the gauge group is 
already unified, these further phases and flavour mixings can succed in contaminating the soft 
breaking matrices before the heavy particle decoupling make them harmless, giving rise to 
important effects [4, 5]. 

Besides the quark and lepton sector and their scalar partner sector, a CP source can come 
from the hidden supersymmetry breaking sector of the theory. As a consequence, there can 
be CP violating effects in the observable sector not associated to a mixing matrix and to a 
consequent flavour violation. As a matter of fact, even in the flavour universality hypothesis, 
a CP violation in the hidden sector can make the soft breaking parameters complex since the 
beginning. In this case, the effects are very large unless the phases are forced to be small. 

In this paper, we consider the effects of phases associated to the misalignment of two Yukawa 
matrices, Au and A D , in the hypothesis that the soft supersymmetry breaking parameters are 
generation universal and real at an high energy scale M . We assume that the Yukawa inter- 
actions are present until Mo without worring about their origin. However the minimal unified 
models that we consider do not give correct unification predictions for the light fermion masses, 
showing that, differently from the MSSM case, flavour physics cannot be totally decoupled 
from unification physics. Furthermore, minimal flavour effects could be comparable to minimal 
unification effects [6]. 

The order of magnitude of EDMs for the quarks u, d and for the electron e depends in a 
crucial way on the gauge structure at the universality scale Mo. We will consider therefore 
different gauge structures, from the minimal one (SU(3)xSU(2)xU(l)) to that one that unifies 
all fermions belonging to the same family (SO(10)), passing through the minimal unification 
(SU(5)), studying them from the point of view of the properties of EDMs under flavour trans- 
formations commuting with the gauge ones. This will allow us to give our estimates for quark 
and electron EDMs. 

We will consider the MSSM in section 2. In Rcfs. [7] the calculation of the down quark 
EDM induced by a loop with charginos has been made using numerical methods with sometimes 
different results. For a given value of tan /3, tan j3 — 10, Bertolini and Vissani claim in fact that 
d n can be as large as d n = O(10~ 30 e cm), Inui et al. obtain d n = O(10 _29 -10 _27 ecm) and Abel 
et al. get d n = O(10~ 33 -10~ 29 e cm), all of them indicating a linear dependence on tan/3. The 
numerical calculation is used in order to solve the matricial renormalization group equations 
for the Yukawa couplings and the soft breaking terms; this is necessary for doing a detailed 
analysis at the Fermi scale. Unfortunatly, in the low energy lagrangian, the particular flavour 
structure of the CKM CP violation, at the basis of the smallncss of the EDMs, is hidden in 
a large number of flavour mixing matrices. On the contrary, by studying the EDMs directly 
in terms of the parameters at the universality scale, it is easy: i) to estimate the order of 
magnitude of the effect, that for tan/3 = 10 does not much exceed 10~ 33 ecm, ii) to display the 
cancellation of contributions not containing the squared light Yukawa couplings, that suppress 
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the EDMs, iii) to show the raise with tan 3 /3 of the d quark EDM. It is also possible to obtain 
the relation dd/rrid = —d s /m s between the one loop dipolcs of quarks d and s, dd and d s (that 
both contribute to d n ). 

From our estimates we see that the one loop supersymmetric effects are comparable with 
the SM ones [8, 9]. It turns out that also supersymmetric contributions at higher loops are 
relevant. 

In section 3 we will consider unified theories. The possibility of EDMs for the neutron and 
the electron close to experimental limits in unified theories has been pointed out and analyzed 
in Refs. [5, 10] in the case of moderate tan/?, where it is associated to the gauge group SO(10). 
The big enhancement of the neutron EDM in comparison with the MSSM case is due to the 
fact that, unlike in the MSSM, in a unified theory CP can be violated also when two eigenvalues 
of the same Yukawa matrix are vanishing [11]. Nevertheless this, as other effects [12], happens 
also in the SU(5) case, but only when the u quark is concerned. As the third generation Yukawa 
coupling involved is the bottom one, the effect is proportional to tan 4 /3. Due to its possible 
interest, we also present the results of a numeric computation of d u . The d quark contribution, 
as the electron EDM, turns out to be very small compared with d u . We conclude section 3 with 
a short rivisitation of SO(10). In appendix, we give the general formulas for the one loop EDMs 
of up and down quarks and charged leptons in terms of the parameters of the Lagrangian at 
the Fermi scale. 

2 EDMs in the MSSM 

In this section we will study quark and electron EDMs in the minimal supersymmetric extension 
of the Standard Model, with real and generation-universal soft terms at some scale M 

A f = Affil, f = u,d,e 
2 2 -, r> /-> 7 j at M . 

m R = m R,0 1 R = Q> u R' d R' L ' e R 

This model provides the minimal amount of CP violation and, consequently, the minimal con- 
tribution to d n in a supersymmetric theory with soft breaking terms generated at M . 

In the described hypothesis, third generation Yukawa couplings cannot be, by themselves, a 
source of quark EDMs. In fact, let us suppose that we neglect the four light Yukawa couplings 
in the quark sector, A M , A c , \d, A s everywhere, a part for the mass that must be present in a loop 
diagram in order to provide the clicity flip. Then, since the hypothesis is scale independent, CP 
violation disappears and the EDMs have to vanish. In fact, at Mp the soft terms are universal 
and real, so that CP violation can only come from Yukawa couplings. But, as there is a couple 
of degenerate eigenvalues in each Yukawa matrix, all CKM phases can be eliminated using 
rotations in the 1-2 sector and phase redefinitions of the fields, just as in the Standard Model. 
More precisely, in the d u case we have to use a rotation of down left light generations, while in 
the dd case we have to use a rotation of up ones. This is because when we consider, for example, 
the EDM of the quark d, the dh mass eigenstate is fixed, and the possibility of rotating the d^ 
and sl eigenstates is lost. In the limit of vanishing A<j and A s , CP violation disappears from 
the theory, but also the lighest down mass eigenstate becomes not defined. 

Thus the light Yukawa couplings play a crucial role in generating the EDMs. Let us now see 
in which way they intervene. To begin with a simpler case, let us consider the imaginary part of 
the B term (defined in appendix) after one loop rescaling from M to Mz- This quantity, like 
the EDMs, vanishes in absence of CP violation and contributes to the EDMs themselves. The 
B term does not depend on the basis in the flavour space in which the left doublets Q, L and 
the right singlets u c , d c , e c are written. In other words, if we consider a transformation of the 
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flavour components of the supcrficlds Q, u c , d c , L, e c commuting with the gauge group, namely 
a U(3) 5 = U(3)q x U(3)„c x U(3)<j<= x U(3)l x U(3) e e transformation, the Yukawa couplings 
become 

\ — > U^c\ v Uq, A d -> UJ c X d Uq, A e — > Uj c X E U L , (1) 

whereas the J5 term remains invariant. For what follows, most important are the U(3)q x 
U(3) u c x U(3)<jc transformations. Since at M all the parameters in the tree level lagrangian 
except the Yukawa couplings are invariant, it is convenient to consider the B term at the Fermi 
scale as a function of those parameters. Actually, the relation between high and low energy 
Yukawa couplings can be inverted, so that the Yukawa couplings can be considered at the Fermi 
scale. 

Owing to the invariance relative to transformations of the right-handed quarks, B depends 
on the Yukawa couplings only through their squares AJA u; a£A d and A E A E . Moreover, if two 
among the A tJ or A D eigenvalues are equal at M , the Yukawa couplings can be made real 
through a transformation (1) so that ImB vanishes. Therefore the RGE corrections to ImB 
must be proportional to (A 2 - A 2 )(A 2 - A 2 )(A 2 - A 2 )(A 2 - A 2 )(A 2 A 2 )(A 2 - A 2 ) « A^A 2 A 2 . 
In fact, owing to the invariance relative to a generic C/(3) 5 transformation, B is a sum of terms 
like 

Tt^XX^KT 1 ■ ■ ■ {\lKr^lK) mr ] Tr[(A E A E ) fe ]< d (2) 

with real adimensional functions of M and A G U d as coefficients. The first non vanishing super- 
symmetric contribution to Im B is then proportional to 

Tr[(A T t A TJ ) 2 (At A D ) 2 (At A L .)(At AJ] (A° u A° d ) ~ A^A 2 A 2 J CP (A° U A° d ), (3) 

where Jqp = ^ m {V dt V tb V^ c V c d) ■ ^ ^-u = A° d , ImB is furthcrly suppressed by lepton Yukawa 
couplings or by small effects due to the iir-cIr hypercharge difference. 

Let us consider now the imaginary part of a flavour non-invariant quantity, more precisely 
the imaginary part of the matrix element of a quantity D transforming as A D (or of a quantity 
U transforming as A tJ ), calculated in correspondence of left and right mass eigenstates. The 
imaginary parts appearing in the expressions (A2) in appendix for the EDMs and thus the EDM 
themselves arc examples of such quantities. In this case, the light Yukawa coupling suppression 
is less strong and dependent on the mass eigenstate we consider, up or down, light or heavy. 

The general dependence of D on the Yukawa couplings is D = A D • /(aJA u; A|,A d , Tr), where 
Tr represents traces of (A E A E ) fc and we suppose that / is a real polynomial in the Yukawa 
couplings, as it is for the one loop EDMs, where the Yukawa couplings only come from vertices 
or supersymmetric RGE corrections. The imaginary part of the d^-dh matrix element is then 
proportional to the d eigenvalue 

lm[D dndL ] = A d Im[/(A T±) dLdL \. (4) 

Moreover, in getting the dependence of f(---)d L d L on the light Yukawa couplings, it is no 
longer possible to consider the limit \ d = A 2 , or A^ = A 2 as above for B, because in this limit 
the eigenstate d,L, as da, is no longer defined. Therefore /(. . . )d L d L has not to be proportional to 
A^A 2 but only to A 2 , so that the necessary dependence on the Yukawa couplings is Im [i^di,] oc 
AjA 2 A 2 Jcp- The first non vanishing contribution to Im[D £ 2 H( i. L ] comes in fact from the term 
proportional to A D (A l t J A u ) 2 (A ] ! ) A D )(AuA u ) in the expression for D. Let a be the proportionality 
coefficient. Since higher order terms are negligible, we have 

lm[D d n df ] ~ a\ dt Im[(At I A u ) 2 (At A D )(At J A u )] dfdf . (5) 
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Had wc considered a matrix U transforming as A U7 we would have found 

Im[U u n uf ] ~ b\ Ul Im[(A ] 5 ) A D ) 2 (A x tA u )(Aj ) A D )] u l u l ■ (6) 

Since the flavour dependence of the one loop EDMs is all contained in the arguments of imagi- 
nary parts in (A2), all the down EDMs are expressable in the previous form = lm\D d rt d L~\ 

with the same D, so as for the up ones we have d Ui = Im[U u R u L ] with the same U. Then there 
arc precise and parameter independent relations between one-loop supersymmetric contribu- 
tions to EDMs of quarks of different families 1 : 

(7a) 
(7b) 
(7c) 

so that 

m _ m i 

(8a) 
(8b) 

Eqs. (8) show that the middle generations EDMs are the largest ones and allow to express 
the neutron EDM in terms of only the u and d quark EDMs. Using eqs. (Al) and (8), and 
neglecting chromoelectric dipole moment contributions, we get in fact the estimate 

d n = xdd + yd Ul with x w 4.6 and y rts 0.7. (9) 

From eqs. (7) it is apparent that dd/md + d s /m s + db/mb is much smaller than dd/rrid- 
This is because the previous combination is invariant under flavour rotations of the quark fields 
and thus suppressed by both X 2 . and X 2 , just as the RGE induced phase of B, electron EDM, 
strong CP angle and 3-gluon operator [13]. Owing to this large suppression, the B term phase 
contributes to the neutron EDM in a negligible way. 

Eqs. (7) show the strong rise of the EDMs with tan/3. The large tan/3 region is therefore by 
far the most interesting one and it is the one that we will consider in the following. Expressing 
all fermion masses in terms of the Yukawa couplings, the weak vacuum expectation value 
v = 174 GeV and tan/3, the only dependence of EDMs on tan/3 is in the Left- Right (and R-L) 
blocks of the scalar mass matrices that provide the "elicity flip" in the scalar sector. In the 
large tan/3 region, 

m drl - - U M D , M 2 RL ~ -v\ v A u (10) 



dd ^ 


+a\d\ t AgA 2 Jcp 


d u ~ — bX u X\X 2 s X^ Jcp 


d s ~ 


— aX s \ t X 2 XIJqp 


d c ~ +6A C A^A^A( Jcp 


db ^ 


+aX b X t X 2 X 2 Jqp 


d t ~ ~bX t X b X 2 s X 2 c Jcp 




d s = d d 

m d 


m b 
= — d b 

m s 




, _ m c 

m u 


m t 
= — d t . 

m c 



so that the EDMs depend on tan/3 only through the down Yukawa couplings: A D oc 1/ cos/3 ~ 
tan/3. Hence the down quark EDMs increase with tan 3 /3, whereas the up quark ones increase 
with tan 6 /3. 

Let us consider now more closely how one loop graphs can give rise to the flavour structure 
of eqs. (7). In both the dd and d u cases there are two relevant graphs (see appendix) involving 
cither charginos or gluinos. The Yukawa matrices necessary to obtain the flavour structure of 
eqs. (5) and (6) come from the RGE corrections to the soft breaking parameters and (in the 

1 We remind that Im(V^bVcV^dV)n = C32[b2i^3i — b3iol2i]Jcp where b,c,d are diagonal flavour matrices, 
b = diag(6i , 62, 63), by = bi — bj , etc. 
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case of chargino exchange only) from the vertices. Each insertion of a couple of Yukawa coming 
from RGEs is accompanied by a loop factor and a large logarithm, tz = (47r)~ 2 ln(Mg /M§). 
In order to provide the 9 Yukawa couplings of eqs. (7), 4 of these inserctions are necessary in 
the gluino diagram, giving a factor t%, whereas 3 are enough in the chargino diagram (t z ). 

To estimate the one loop contributions to dd and d u only the dependence on dimensionful 
parameters is missing. Because of the necessary presence of a L-R scalar mass inserction as in 
cq. (10) and of the behavior of the loop functions, it is 

4\2 X 2 v ^ A u T _ 1fl -3i Aan/3\ 3 /200GeV\ 2 



~ h > V 55 / V M SUSY / ( Ua ) 



4 « J^tlXdXtXtXl ( Z A \, 2 Jcp « 10" 31 ecm i 

(4tt) 2 max(m^,m^)m^ 

"rf ~ e -r t zXd,X t X b X c — ~ ST — 5- Jqp ~ 10 ecm 

d 4 z c max(M|,m?)m 2 ; \ 



3 /onn n ia2 



"3 +4 x ,4 ^2 X 2 ^3 7 __ _ , n _ 32 „ _ Aan/3\ ( 20 ° GeV 

d>" l d v 55 / V^susy/ (nb) 

for the chargino and gluino one loop contribution to dd respectively. In this and in the following 
estimates, we neglect all the numerical coefficients of order one, we choose M « 10 16 GeV and 
we use central values for the various known parameters, in particular, Jqp — 2 • 10~ 5 . Moreover 
'tosusy' stands, in each case, for the particular combination of soft parameters written in 
the analytical approximation. Because of d u /dd ^ (tan/3/55) 3 and of x/y w 7 in eq. (9), the 
corresponding contributions to d u arc less interesting and totally negligible for moderate tan (3. 
At a closer inspection, these estimates turn out to be correct within an order of magnitude. 

As it is apparent from eq. (11), d n cannot be much larger than 10 _33 e cm for tan (3 = 10 and 
than 10~ 31 ecm for whatever value of tan/3. We remark that, owing to the necessary presence 
of A 2 in all results for dd, a numerical calculation of the d quark EDM requires the knowledge 
of the low energy parameters with a precision of about 1/10 6 in order to see the reciprocal 
cancellation of the terms not proportional to A 2 . 

It is at this point interesting to understand why the three-loop pure SM contributions to the 
neutron EDM, rf^ M sa 10 _32 ecm [8], are not suppressed relative to the one-loop supersymmetric 
ones. 

A large suppression of the pure SM contribution would seem plausible because, in comparison 
with it, the supersymmetric one loop contributions can be enhanced by 

1. large logarithms, log Mq/M|, that multiplicate the loop factors (47r)~ 2 in the RGE in- 
serctions involved by the one loop supersymmetric diagrams; 

2. tan/3 factors, that can enhance the A D and A E couplings; 

3. a particularly small effective combination of SUSY parameters, 'tosusy'- 

On the other hand, the SM contributions are enhanced in different ways. In fact, in the 
pure SM graphs, 

i. the elicity flip can occur on an external leg, giving a factor equal to the "constituent quark 
mass" m™ nst - m„/3, equal for q = {u, d, s} (so that df* w <if M ); 

ii. the Yukawa couplings are expressed in terms of quark masses evaluated at low energy and 
therefore enhanced by QCD renormalization factors, as the strong coupling; 

iii. unlike the RGE-induced contributions, that must depend in a polinomial way on the 
Yukawa couplings, the infrared structure due to subtractions of quark propagators in the 
SM graphs, give rise to a mild logarithmic dependence. 
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Point iii. is particularly important in the case of quantities neutral under the quark U(3) 3 
flavour group (like the electron EDM, or the phase of the i3-term, or the strong CP angle, or 
the 3-gluon operator [13]) that for this reason receive purely supersymmetric contributions much 
smaller than SM ones. Due to the absolute lack of experimental interest, we avoid discussing 
such issues in any detail. 

Other than one loop supersymmetric and three loop pure SM diagrams, it is also possible 
to have higher loop contributions, not considered before, having the two kind of enhancement. 
An example can be obtained from the chargino one loop diagram by adding by adding a QCD 
loop, that allows elicity flip on the external lines and gives a contribution to dd comparable 
with the previously computed one. Using charged Higgs exchange, it is also possible to obtain 
an interesting three loop contribution to d n given by 

, H ± a 3 iTLm?,,,,, „ , 9 /tan /A 4 /200 GeV\ 4 

(4tt) 5 m% ± \ 55 J V m H± J 

Fine tuning considerations suggest that, in the large tan/3 region, the charged Higgs are lighter 
than the squarks [14]. 

To conclude this section, let us note that the universality scale Mq can be much smaller 
than the unification scale if the supersymmetry breaking is transmitted at relatively low energy 
by gauge interactions. In this case all RGE-induced contributions to d n are smaller because of 
their strong dependence on tz 1 whereas the higher loop ones, and in particular (12), are not. 

3 EDMs in unified theories 

In this section, we suppose that the gauge group is unified at the universality scale Mo that 
we will identify with the reduced Planck mass. In this case, EDMs can be much larger than 
in the MSSM. As seen in the previous section, the smallness of EDMs in the MSSM is due to 
the light Yukawa couplings, that have to be present to prevent the removal of all phases from 
the lagrangian, phases that at the Planck scale reside only in the Yukawa matrices. In fact, 
with regard for example to the quark d, if \ u = \ c = CP violation can be removed with 
independent redefinitions of the right down quark and of the left doublet. If the gauge group is 
unified, at the Planck scale the phases still reside only in the Yukawa matrices, but in this case 
MSSM multiplets belonging to the same representation of the unified group can no longer be 
rotated independently. Hence, depending on the gauge group and on the quark in consideration, 
CP violation can persist even if some or all of light Yukawa couplings are vanishing, making 
not necessary their suppressing presence in the EDM expressions. 

Let us consider first minimal SU(5) unification. One generation is composed by two multi- 
plets, a five-plet F — (d c , L) and a ten-plet T — (u c , e c , Q), and the Higgs fields belong to two 
five-plets, H, that transforms like F, and H, that transforms in the conjugate way. In terms of 
these fields, the Yukawa interactions are 

^TiX Ui TiH, V^FiX^TjH (13) 

in a basis in which X XJ is diagonal. Below the unification scale, A (J and A D are the same of the 
previous section. 

Five of the six phases of the CKM matrix can be removed in the MSSM through inde- 
pendent redefinitions of the Qi, u^, d\. Here, it is easy to show that only three phases can 
be removed by independent transformations of Tj, Fj, so there are two more physical phases 
compared to the MSSM (and the SM). Whereas two phases decouple below the unification 
scale in non supersymmetric models, in supersymmetric models they leave their effects in soft 
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breaking parameters, that loose their universality and reality, before decoupling could make 
them harmless. 

Is the effect of these two further phases on the EDMs suppressed by light Yukawa couplings 
in this model? The answer is no for the quark u and yes, but in a different way compared to the 
MSSM, for down quarks and charged leptons. In order to show that, let us suppose first that 
the two light A„ eigenvalues, X u and A c , are vanishing. In this case, all phases in the Yukawa 
matrices become unphysical. In fact, the first of the two interactions (13) becomes T 3 X t T 3 H /4 
and it is invariant for transformations in the T x -T 2 sector. Moreover, as it is well known, it is 
possible to write the most general unitary matrix V* as 



T/t = diag(e»,e^,e-)(-^^) 











R 2 , 



R 3 











diag(e M ',e l ' 3 ',l), 

(14) 

where R\, R 2 , R3 are orthogonal 2x2 matrices. Then we can reabsorb e M ,e l ^,e l7 ,e M , e 1 ^ 
by a redefinition of F 1 ,F 2 ,F 3 ,T 1 ,T 2 , reach the phase e iS by defining (T{,T^) T = R 3 {T 1 ,T 2 ) T 
and reabsorb also e lS by a redefinition of T[, without affecting the diagonal interaction. Since, 
as seen in the previous section, the limit X u = X c = is meaningful for the down quarks and 
for the charged leptons (and for the top), this means that their EDMs are suppressed by light 
up Yukawa couplings in some way. 

Let us suppose now that the two light A D eigenvalues, Xd and A s , are vanishing. In this case, 
two phases remain physical. In fact, the interactions (13) become 

i TiX Ui Ti H, V2 FsXbV^Tj H. 

With a redefinition of F3, one of the three phases of V 3 \, V 32 , V 33 can be reabsorbed, but not 
the remaining two because, if they were reabsorbed by a redefinition of some Ti, they would 
reappear in the diagonal interaction. The two remaining phases are just those ones related with 
unification. Since the limit A^ = X s = is meaningful for up quarks (and the bottom and the 
tau), their EDMs are not necessarily suppressed by light Yukawa couplings. 

The dependence of EDMs on Yukawa couplings can be obtained as in the MSSM. Changing 

basis in the flavour space for the SU(5) supermultiplets T, F corresponds to make a U(3)t x 
U(3)-p transformation on the flavour component of the fields, relative to which the Yukawa 
couplings transform in this way: 

A LI - U t T KUt, A d - U ¥ T \ D U T . (15) 

As before, the down quark EDMs are given by dd t = ^[-^d^df ] 1 where the matrix D 
transforms as A D . Unlike the case of the MSSM, this does not mean that D = A D /(A| J A U , A£A d ). 
In this case we have rather 

^ = A D /((AtA D ),(AtA D )*,A u ,A* I ), (16) 

where / depends on the arguments in such a way that / — > U^fU T for a transformation (15). 
From eq. (16) it follows that dd is proportional to A^ as in the MSSM. Moreover, / must depend 
on down Yukawa couplings through their squares but it can depend on individual up Yukawa 
couplings also through not squared couplings, provided that the total number of A lT is even. 
This can be seen by noting that the particular flavour transformation T — ► iT, F — > F leaves 
AnA D and / unchanged but it changes sign to A TI . Therefore, unlike the MSSM, a A c suppression 
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can be enough for the quark d. Indeed, the first non vanishing contribution to dd is proportional 
to 



Im[(A D At(AtA D )*A^) dRd J = X d X 3 t X 2 X c Im[vlV cb v£v td \ . (17) 

The EDMs of the quark s and of the electron are suppressed in an analogous way. 

The up quark EDMs behave in a different way. They are given by d Ui = Ivci\U u r u l\ where 

the matrix U transforms as X v under the transformation (15). Unlike the MSSM, U is not 
necessarily proportional to A tJ on the left, but can also be proportional to A^ , that transforms 
in the same way on the left. Therefore, d u is not necessarily proportional to A u; giving rise to 
a possible enhancement. In any case, d u must contain an odd number of up Yukawa matrices, 
since the transformation T — > iT, F — > F leaves A D unchanged but it changes sign to A LI and 
U. The first non vanishing contribution to d u is proportional to 

Im[(At \ D )*\ v (\i A D )] = XX lm[V uh V b T t V tb V b {] (18) 

and it exhibits the described features. In eqs. (17,18) we omitted the renormalization factors 
for the Vckm matrix elements. 

In SU(5) we have therefore d u 3> dd,d s ,d e . Only the gluino diagram is able to give rise 
to the At enhancement. In fact, the chargino diagram is explicitely proportional to the X u 
that appears in one of its vertices and it is negligible in this model. Let us concentrate then 
on the gluino contribution to the quark u, d\. The flavour structure (18) is generated by the 
corrections to the mass matrices. At the unification scale, these corrections are proportional 
to combinations of Yukawa matrices transforming as the mass matrices themselves relative to 
the U(3)t x U(3)-p group. Hence the corrections to the ten-plet mass matrix, and so those 
to rrtQ, rn 2 * R , m 2 * R , are proportional to 1, AJ,A u; a£A d , etc., whereas the corrections to the 
five-plet mass matrix, and so those to m\ and rn 2 * da , are proportional to 1, (A D An)*, etc.. All 
these corrections are also proportional to t G = (Air)~ 2 \og(Mp / M G ) ~ 0.06. Some of them 
are characteristic of unification, whereas other ones are produced also below the unification 
scale and hence at the Fermi scale they are proportional to tz — (4-7r) -2 log(Mp/M§) ~ 0.5. 
Including also some approximate numerical factors, the corrections Am 2 at the Fermi scale are 

Am| oc 1 - 3t z Xl\ i; - 3t z XiX D (19a) 
Aml R ocl- 6t z X v Xl - 6t G (XlX D )* (19b) 
Ato 2 ^ oc 1 6t z X D Xl. (19c) 

These corrections are able to generate the flavour dependence of eq. (18), so that we can estimate 
\d n \ r, e—Xtlm[V ub V bt V tb V b l\tzt G 



in-(25-^26) 

/tan/3\ 4 f \Jm\V%^ \ { 500 GeV 
I 50 M 10- 5 J V^susy. 



(20) 

e cm, 

where Msusy in the second line is the combination of supersymmetric parameters appearing in 
the first line and it can be even larger than 500 GeV because of the likely lightness of A in the 
large tan (3 regime. Since the prediction for d u is interesting enough, we use the expressions (A2) 
to do an exact computation of d n that we show in figure 1 for M 3 = 500 GeV, for three values 
of tan/3, 2, 10 and 50, and for | Im[^F ub ]| = 10~ 5 . d u is plotted as a function of M3 / m? in 
order to exhibit the so-called "gluino focussing" effect [10]. 
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Figure 1: The contribution to the neutron EDM generated by minimal SU(5) for M3 = 500 GeV ? 
tan j3 = 2, 10, 50 and | Im[V^|V u( ,] | = 10~ 5 , as a function of the ratio between gaugino and scalar 
masses and for random samples of acceptable sparticle spectra. 



From eq. (20) and Fig. 1 and from the experimental limit d n < 0.8 • 10~ 25 ecm [18] it is 
apparent that the u quark EDM can be very large in superunified theories with large tan (3. 

Let us consider now shortly the gauge group SO(10), whose sixteen dimensional representa- 
tion "16" unifies all quarks and leptons of one generation, included a right-handed neutrino. In 
the minimal model the two light Higgs doublets belong to two different ten dimensional SO(10) 
multiplets, 10„ and 10^, and the Yukawa interactions are 

16 T A tJ 16 10„, 16 T A D 16 10 d . (21) 

In this model the situation is different from the SU(5) case, in which the five-plets and ten- 
plets could be rotated independently. All six phases of V are physically significant. Moreover, 
it is easy to see that in this case it is not possible to remove the CKM phases neither when 
Ad = A s = (like in SU(5)) nor when X u = \ c = (unlike in SU(5)). Therefore neither d u nor 
dd are suppressed by light Yukawa couplings. 

The u quark EDM is generated in the same way it is generated in SU(5) and its estimate is 
identical. On the other hand, whereas in SU(5) only mf lR (and not m d ) has non MSSM cor- 
rections, in SO(10) the situation is symmetric under u-d exchange, so that rn dR has corrections 
depending on A D as the rnf lR ones depend on X v . Moreover, also an approximate expression 
for dd can be obtained from that one for d u exchanging u and d. dd has therefore a linear 
dependence on tan/3. Moreover, when tan/3 is large and Xb — At the main differences between 
d u and dd are due to the different combination of CKM angles and phases, Im[Vj d V^;] instead 

of Im[V^ 6 V^], and to the different RL flip, that is provided by /i and not A u . With similar 
phases, the CKM angles could favour dd- Nevertheless, also in SO (10) with large tan/3 the u 
quark contribution cannot be forgotten. 

It is also interesting to estimate the corrections to the quark masses that, in SU(5) and 
in SO(10), contribute to the strong CP violating angle 9 in models without an axion. These 
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MSSM 


minimal SU(5) 


minimal SO(10) 


d u 
dd 

d e 


A„A^A( Jcp 
A d A* A^A^ Jcp 
•^eAj A^Ag A^ Jqp 


\ t xlMv 2 ub v t l] 
XdXlx 2 b XcMv cd v cb v tb v t d} 

X e X 3 t XlX c lm[V cd V cb V tb V td } 


X t Xilm[V Z ub V t l} 

x b xtMV 2 td v t l] 

X T Xflm[V 2 td V t l] 



Table 1: Flavour factors suppressing the purely super symmetric contributions to light fermion 
EDMs in different minimal models. 

contributions are 

^cd - ^tHVl b V t l]t z t G « 10--tan^(«^H) (22) 
both in SU(5) and SO(10) and 

«U - £^Hfett*$ - ,23) 

in SO(10). The bound 0qcd < 10~ 9 shows that an axion is necessary in SO(10) and also in 
SU(5) if tan/? is large. It is interesting that the large amount of CP violation left by SO(10) 
unification in the soft terms furnishes an experimental possibility to see the 'invisible' axion [19]. 

4 Conclusions 

We have studied EDMs produced by the misalignment of two Yukawa matrices in models 
with universality of soft breaking terms at an high energy scale from the point of view of 
transformation properties under flavour transformation commuting with the gauge group. In 
this way it is simple to get the dependence of EDMs from Yukawa couplings, summarized in 
table 1. On this basis, estimates are possible and effective for small and large effects. 

In the case of MSSM, owing to light Yukawa coupling suppression, the effects given by 
one loop diagrams are largely below present experimental limits (see eq. (11))- Therefore, the 
detection of a non vanishing EDM for the neutron in foreseeable experiments would be an 
inequivocable signal of physics beyond the MSSM. Some interesting contribution to d n can 
arise also from more loop diagrams. 

In the case of minimal supersymmetric SU(5) and large tan (3, on the contrary, the effects 
can be close to experimental limits with regards to neutron EDM because of the contribution 
of u quark, proportional to tan 4 /3. For this interesting case the results of an exact computation 
are given by Fig 1. The smallness of d quark and electron EDMs compared with the u quark 
one is also characteristic of SU(5). Therefore SU(5) can be distinguished from the MSSM 
because it can give rise to a measurable EDM for the neutron. Moreover, SU(5) effects can be 
distinguished from the SO(10) ones and from effects coming from universal intrinsic phases in 
A and B terms [20] because of the smallness of the electron EDM. On the other hand, at least 
for the EDMs a similar pattern can be given by strong CP violation. 

In the case of minimal supersymmetric SO(10) it is well known that the d quark and electron 
EDMs can be close to the experimental limits. If tan/3 is large, the effect is enhanced by a 
factor /ita,nf3/Ad- Also in this case, as in SU(5), the u quark gives an important contribution 
that must be taken in account. 
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Appendix: general formulas for one loop EDMs 



It is possible to express the neutron EDM d n in terms of the contribution of the quark q to the 
neutron spin (Aq) n and of its EDM d q [21]: 

d n = ((Ad) n d d + ((Au) n d u + ((As) n d s , (Al) 

where (Ad)„ = 0.82 ± 0.03, (Au) n = -0.44 ± 0.03 e (As) n = -0.11 ± 0.03 [21] and C is a 
renormalization factor, ( ~ 1.6. 

Concerning the quark EDMs, let us consider a supersymmetric extension of Standard Model 
at the Fermi scale, with minimal field contents, imposed R-parity and generic soft terms a part 
from a reality hypothesis on B (defined as the coefficient that appears in the soft breaking term 
B/j,H u Hd, where (i is the coefficient of the corresponding interaction in the superpotential) from 
which it follows that the chargino's and neutralino's mixing matrices can be taken real. Such 
hypothesis is verified with very good precision in the models that we will consider. Then the 
full expressions for the one loop EDMs of up and down quarks are 



d Ui 



(47T) 
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8 a 3 1 
+ 3 6 4^M3 
v(ur) 



nh d M x 
Im ' 
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* 3 (dL) 
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A„ 



M 2 

Id )A\ 



where £3 (a) and y(a) are respectively the third component of weak isospin and the hypcrchargc 
of the particle a (normalized in such a way that q = t 3 + y). g is the gluino and M 3 its mass, 
N n , n = 1 . . .4 are the neutralinos, x«> Xn > n = 1 • • - 2 the charginos and H~ their 

mixing matrices. M^, M|, are the 6x6 squarks mass matrices in the mass eigenstate basis of 
corresponding quarks 



M 2 = 



Ml = 



' 1 ( u 
mi 



Du L l 



M 2 

M U (A U + ft cot 



"I 

-M d (^ a 



M 2 



fj, tan /? 1) 
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-(At + MCot/31)M„ N 
-(A^ + ^tan/31)M d N 



(A3) 



In the right sides of (A2) we used a matrix notation. The fs are the appropriate loop functions, 
namely linear combinations q s g2 where q s and qj are the electric charges of the particles 

running respectively in the scalar and fermion line of the corresponding diagram (in unity of 
e), and 

= 01 1 lVa [r 2 -2rlogr-l], h 2 {r) = - 1 [-2r 2 logr + 3r 2 - 4r + 1]. 
z \ r - L ) z \ r - L ) (A4) 

In eqs. (A2) the first two contributions to EDMs are due to charged higgsinos and gluinos 
exchange and they are the dominant ones. The third ones take account of bino and neutral 
wino exchange and they are less important than the corresponding gluino exchange, while the 
last ones come from neutral higgsino exchanges and they are completely negligible. All the 
contributions in eqs. (A2) come from diagrams with none or two elicity flips on the vertices. 
The one loop diagrams in which the elicity flip occurs on an external leg are always real. As 
such they do not give any contribution. The diagrams with one elicity flip on the vertices would 
contribute only if the B term were complex. In this case in eqs. (A2) also the mixing matrices 
H would be complex and they would appear within the imaginary parts. 

The arguments of imaginary parts contain both the model dependence and the difficulty of 
the calculation of EDMs. 
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